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Abstract: 



Based on a sumrule for the nucleon spin we expand quark and gluon or- 
bital angular momentum operators and derive an evolution matrix for higher 
moments of the corresponding distributions. In combination with the spin- 
dependent DGLAP-matrix we find a complete set of spin and orbital angular 
momentum evolution equations. 
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1 Introduction 

Since several years much effort is invested into the determination of the in- 
ternal nucleon spin structure by means of polarized lepton-nucleon deep- 
inelastic-scattering. One of the focal questions in this respect is the deter- 
mination of the gluon spin fraction AG(x, Q 2 ), which could be as large as 2, 
though its actual value is only poorly known. Let us refere here just to some 
recent analyses [0, [| f|, ||, § instead of reviewing the slightly involved history 
of this discussion. As the total angular momentum of the nucleon is one half 
and the quark spin is known experimentally to be small, a large gluon spin 
implies substantial angular momentum. It is a natural hope that angular 
distributions of final state hadrons can be linked to the angular momenta of 
quarks and gluons inside of the struck nucleon. Up to know it was, however, 
not possible to establish this connection within a self-consistent QCD scheme. 
This would require e.g. the correct treatment of transverse momentum, in- 
cluding all spin-momentum-correlations, during the fragmentation process. 
Here we want to disregard such problems and focus instead on the com- 
bined evolution equations for quark and gluon spin and angular momentum 
distributions. This evolution was studied for the first moment quite exten- 
sively by e.g. Xiangdong Ji We extend his analyses to higher moments. 
The relevance of such investigations for the possibilities to find experimen- 
tal observables linked to the partonic angular momentum is probably rather 
obvious. The x and Q 2 dependence of the angular momentum distributions 
will be a crucial ingredient to decide whether it is possible to find for ex- 
isting or planned experiments an x-Q 2 -window for which simultaneously Q 2 
is large enough for perturbative QCD to apply, v is large enough to assume 
leading quark fragmentation, and the spin-asymmetries are so large, that an 
experimental determination seems possible. 

2 Preliminaries 

Let us begin with a presentation of a partition formula for the spin of a nu- 
cleon in an helicity-eigenstate in light-cone coordinates (LCC). The angular 
momentum density is defined through 



which is a gauge invariant rank-3 tensor. Let now \P,S) be a nucleon-state 
with the normalization 





(P', S'\P, S) = 2P+(2vr) 3 5 3 (P - P')6ss>- 



(2) 
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In the following we consider the quantity 

M^ X (P,P',S) = J d 3 x{P',S\M fU/X {x)\P,S) (3) 

and treat it as a distribution in the variable P' (in LCC, the integral becomes 
/ dx~ dx 1 dx 2 ) . To avoid derivatives of delta distributions or delta distribu- 
tions at the point 0, a bilinear functional 

(7(P, P'), $(P')) = / d 3 P'f(P, P'MP') (4) 

is defined with a suitable test function $(P') for an arbitrary distribution 
f(P, P'). A calculation quite similar to the one in |12[ leads to 



(5) 



M+ 12 (P,P',S),<5>(P')) 1 



2P+(27r) 3 $(P) 2' 
where = (P+, 0, 0, P~), ^ = (P+, 0, 0, -P~), and 

'(P', ^|P, S), $(P')) = 2P+(2tt) 3 $(P) (6) 



plays the role of the norm (0) . In the case of QCD, the above component 
of the angular momentum density tensor is given by (in the following we use 
the light cone gauge A + = 0) 

M +l2 (x) = i^ 1 + (x l d 2 -x 2 d 1 )^ 

+d + A j (x 1 d 2 -x 2 d 1 )A j 

+A l d + A 2 - A 2 d + A\ (7) 

where j is summed over 1 and 2; and d + = j^- = = cL. Using the 
definitions (see Jl|, |I2]| ) 

AS = {P\S\Jd 3 x~* 1 +[ 1 \ 1 2 ]*\P,S) (8) 

Ag = (P' 1 S\Jd 3 x(A 1 d + A 2 -A 2 d + A 1 )\P 1 S) (9) 

L q = (P',S\ J d 3 xiU 1 + {x 1 d 2 -x 2 d 1 )V\P,S) (10) 

L g = (P',S\Jd 3 xd + A j (x 1 d 2 -x 2 d 1 )A 1 \P,S) (11) 



and (|5]), one gets 

(AZ/2 + Ag + L q + L g ,®(P')) , 

= (12) 

2P+(2tt) 3 $(P) 2 1 ; 

This is the nucleon spin sumrule. Because the 3-space integral over M +12 (x) 
is equal to the third component of the generator of rotations, (|T2]) can be 
seen as a decomposition of the nucleon spin into quark and gluon spin and 
orbital angular momentum contributions [I0|. Ag turns out to be the first 
moment of the difference of the gluon distributions with longitudinal spin 
projection ±1. 

Ag(t) = dx\g^(x,t) - gfat)]. (13) 



Due to the general scale dependence of the theory, the quark and gluon 
distributions are functions of t = ln(Q 2 / Aq CD ) (Q 2 being the typical scale 
set by the external momenta). The t evolution of the n-th moments of the 
spin dependent parton distributions, Ag n (t) and AE n (t), is given by the well 
known DGLAP-equation (0,0). The evolution matrix is the anomalous 
dimension matrix of the corresponding operators, e.g. for the quark spin 
operator (in normal coordinates) 

S q = J d^^yV,-^]* = / d 3 x^ 7 3 75^ 
— , §»i...»n-i = J d z x^^^D^ . . . r^- 1 ^, (14) 

which are the familiar operators of the operator product expansion for deep 



inelastic scattering (see for example JT4 



In view of our aim to calculate the t evolution of the functions L q {t) and 
L g (t) generalized to higher moments, L™{t) and Lg(t), we expand by analogy 
the quark and gluon orbital angular momentum operators: 

L q = J d 3 xi^- 1 + {x 1 d 2 -x 2 d 1 )^> 

— ► £w-Mn-i = J d 3 x i^ 1 + (x 1 d 2 - x 2 d v )D^ . . . D^- 1 ^ (15) 

L g = J d 3 x d + A ha (x 1 d 2 - x 2 d r )A j ' a 

— > ^1-/^-1 = J d 3 xd + A j ' ai (x 1 d 2 - x 2 d 1 )D^ 1 ' aia2 . . . D^- 1 ' an - ian A j > an 

(16) 

is the covariant derivative <9 M — igT a A^' a (respectively D^ ab = 5 ab d^ — 
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3 Calculation of the evolution coefficients 



Due to operator mixing, we have to calculate the 4x4 anomalous dimension 

matrix for the set of operators J = (S 1 ^ 1 "^" 
which is given by 



Si* 1 - 

' ^9 



•Mn-l r^l.-./in-l T [il-..fJ,n 



la 



7-1 



dZ 



d(Zol 



<91n(A 2 



d In A 2 



+ 0(g 4 



(17) 



where the matrix Z denotes the operator renormalization matrix, defined 
through 



- ren - unren 

J t = (Z~%J J 



To determine Z 1 we inspect the insertion of our generalized orbital angular 
momentum operators (|Tj| ^) into off-forward matrix elements 



(19) 



(\p, s) is a quark or gluon helicity state) and calculate their divergent parts. 
In the following, we restrict ourselves to LO, i.e. only 1-loop-graphs are con- 
sidered, and work entirely in LCC and the LC-gauge In an interacting 
theory, the momentum and spin eigenstates can be expanded pertubatively 
order by order in the coupling constant g 



\p,s) 



\P, s) + 
\p,s)o 



-H + ie 
\p,s)i + 



H int \p,s) + O(g 2 ) 
0(g 2 ). 



(20) 



Here is just the quark- gluon or three-gluon interaction vertex. As we 
are only interested in the twist-2 part of the Q 2 -evolution we have only to 
deal with the (/ii = +,..., // n _i = +) components of the operators. Graphs 
like figure le, which come from the covariant derivatives in (|1|,[T^), do not 
contribute in the light cone gauge []16|. An explicit calculation shows that 



the matrix elements have the following general structure (here F(p,p'; n) and 
G(p,p';n) just stand symbolically for the individual contributions listed in 
Eq. (26)-(29)) 



(p',s\L+ g -- + \p } s) ll <S>(p') 



(27T) 



d d 
ip 2 ^ i -ip 1 ^)5\p'-p) 



dp 1 '- 



F(p,p';n),$(p') 
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and 



= " ((2tt) W - p) (V ^ - ip^y&p'; n), d>(p')) 

- ^27rf5 3 (p> -$)F(p,p';n), (ip 2 -^ - ip 1 ^)*(po) • ( 21 ) 
Observing that 

'(p / , S |L+ 9 + b, S ) ,$b / ) / 
- ((-^ + )"- 1 2p + (27T) 3 ( 5 3 (p"-p), (y^r-zp 1 ^)*(P')) ( 22 ) 

((2vr) 3 5 3 (p' - p) (ip 2 ^ - i p 1 ^)F(p,p>; n), $(p')) = 

= (G(p,p»(p', S |5+ 9 - + b, S ) ,$b')) (23) 

the elements of Z~ x are easy to determine: The first term on the r.h.s. of 
Eq. fl2T|) gives the contribution from the spin operators and the second the 
contribution from the orbital angular momentum operators to the renormal- 
ization of L+ g ' + . 
With the definitions 

(Op) 9 o 9 = ((p',s\Op\p,s)r ark ' 9luon ^(p')) (24) 
(Op)™ = ((p\smp,s)r rk ' 9luon ,Hp')) (25) 
the divergent parts of the matrixelements are (figures la-Id) 



CV / d:rx n (- ) — Cf I dx(- 



o 1 — x Jo 1 — x 



o 



+ C F [ 1 dxx n - 1 (x 2 -l){S+-+) q \, (26) 
^ 2 



(L+-+)? = _ln(_){ n/ jf dxx"-V + (l-^)(^r + )o 

+ n / ^ 1 rfxx n - 1 (-2x 3 + 4x 2 -3x + l)(5+"' + )o}, (27) 

(^■■ + )? = ^^(^{^/^^^(l + a-x) 2 )^)^ 

+ C F J dxx n -\-x 2 + 3x-2)(5+- + )g|, (28) 



a A\ 

Z7T /T 



2C A / dxx r 



_ 1 x 4 -2x 3 + 3x 2 -2x + l 



1 — x 



^ f 1 dx (2x 2 - 2x + 1) 
2 jo 



9a / 1 da;( ^±^_ (a; 2_ x _ 1)) 

2 Jo V 1 — ar V ;; 



(Lf + )t 



+ 2CUjf dxx n -\x- \){x 2 -x + 2)(S+-+)l 



(29) 



A is the UV-cutoff from the transverse momentum loop integral (li is the 
infrared cutoff). The terms a, b and c correspond to the quark respectively 
gluon self energy, coming from the corresponding field renormalization con- 
stants Z~g which are implicitely included in the matrix elements. In combi- 
nation with (|T7|) and the definition ^-Ay = 7y, we obtain the generalized 
and completed evolution equation for the n-th moments of spin and orbital 
angular momentum distributions 



d 
dl 



T " 



A 



ss 



A 



SL 



A 



LS 



a(t) ( A 
2tt [ A 



ss 

n 
LS 



An 
An 



\L n Q J 



O sr^n 1 
Z ^.7=1 n 




( AS" ^ 



T " 



\L n Q J 



2C, 



n +n+6 



n .f 



n-1 



n(n+l) 
11 _ rif ' 2 _ 
1"? 6C A n(n+l) 



n(n+l)(n+2)(n+3) 
n 2 +4n+6 



n(n+l)(n+2)(n+3) 



■3=1 



2C, 



n 3 +3n 2 — 6 



n 2 +3n+4 



n(n+l)(n+2)(n+3) 



(n+l)(n+2)(n+3) 
QC A ^ 
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is the known DGLAP-matrix of the spin dependent case. To our knowl- 
edge the coefficient matrices A^ s and A\ L were not calculated before except 
for n = 1. A^ L turns out to be equal to the spin-independent DGLAP-matrix 
withn->n + l (@,[[7|). 



The results for Ag L and A\ L can be interpreted in the following manner. 



Angular momentum can only be sensibly defined between at least two par- 
ticles. On the level of the twist- 2 DGLAP equations which describe the 
splitting of one parton into two the angular momentum of the initial parton 
cannot play any special role. It is only relevant for the final state, which 
contains two partons. Here it allows for a spin flip. The momentum (longi- 
tudinal and transverse part) of the incoming parton is just split between the 
daughter partons with the ratio x/(l — x). Therefore A\ L turns out to be 
the usual spin-independent results just shifted by 1 in n, because the angular 
momentum operator contains the momentum as a factor. 
Following the same arguments one has to conclude that the angular momen- 
tum carried by the initial parton (relative to the nucleon center of momen- 
tum) does not couple to the final state spins. Therefore A$ L has to be zero. 

To check our results we repeated our calculation following the manner in 
which the splitting functions P(x) are usually derived and just inserting the 
angular momentum operators Eq. flSD-fllll). In this case some care is needed 
when the singularities of the diagonal splitting functions are replaced by 
the (yz^)+ distributions 0. Let us note that Eq. ( |T2"D implies 

d (AE(t)/2 + Ag(t) + L q (t) + L g (t), $(P')) 

_A L — n C\\\ 

dt 2P+(2tt) 3 $(P) 1 ; 

which is just another check. 

The estimation of the phenomenological implications of our results requires 
numerical simulations with a great variety of initial angular momentum dis- 
tributions. These studies will necessarily invoke a substantial amount of 
model-building, while the results presented in this contribution are valid in 
general. 
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